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Abstract 

Bipolar disorder is characterized by internally affective fluctuations. The abnormality of inherently mental state can be 
assessed using resting-state fMRI data without producing tasl<-induced biases. In this study, we hypothesized that the 
resting-state connectivity related to the frontal, striatal, and thalamic regions, which were associated with mood regulations 
and cognitive functions, can be altered for bipolar disorder. We used the Pearson's correlation coefficients to estimate 
functional connectivity followed by the hierarchical modular analysis to categorize the resting-state functional regions of 
interest (ROIs). The selected functional connectivities associated with the striatal-thalamic circuit and default mode network 
(DMN) were compared between bipolar patients and healthy controls. Significantly decreased connectivity in the striatal- 
thalamic circuit and between the striatal regions and the middle and posterior cingulate cortex was observed in the bipolar 
patients. We also observed that the bipolar patients exhibited significantly increased connectivity between the thalamic 
regions and the parahippocampus. No significant changes of connectivity related to the frontal regions in the DMN were 
observed. The changed resting-state connectivity related to the striatal-thalamic circuit might be an inherent basis for the 
altered emotional and cognitive processing in the bipolar patients. 
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introduction 

Bipolar disorder is a mood disorder exhibiting a prevalence of at 
least 1%, constituting a considerable health care burden [1,2]. 
Internally affective disturbances of the mental state, including 
mania and depression, are the principal symptoms of bipolar 
patients. Patients in a manic state typically experience an increase 
in energy, racing thoughts, and a decreased need for sleep. In a 
depressive state, patients experience sadness, guUt, hopelessness, 
disturbances to sleep and appetite, and a loss of interest in usually 
enjoyable activities. Additionally, bipolar patients have been 
reported to exhibit deficits in several cognitive functions, such as 
sustained attention, executive function, verbal memory, and 
decision making [3-6]. 

Emotional and cognitive processing, which are impaired in 
bipolar patients, are mediated by the frontal-striatal-thalamic 
(FST) circuit [7-9] . The FST circuit involves in diverse functions, 
including mood regulation, reward processing, action selections, 
strategic planning, and working memory. The functions of the 
FST circuit rely on the various interactions among the regions in 



the circuit [8,10, 1 1]. For example, the thalamus participates in the 
FST circuit for relaying striatal inputs to the frontal regions and 
providing feedback to the striatum [7,8]. Additionally, the medial 
frontal regions of the FST circuit are connected with the posterior 
cingulate cortex and the temporal regions in the resting state and 
these regions resemble the default mode network (DMN) which is 
implicated in emotional and self-referential processing [12]. The 
interactions between the FST circuit and the DMN regions may 
also be involved in the pathophysiology of bipolar disorder [9] . 
Although previous neuroimaging studies have reported structural 
and functional abnormalities in the regions in the FST circuit of 
bipolar patients [13-17], the interaction changes within the FST 
circuit and between this circuit and other regions have been 
explored less thoroughly. 

Interaction between brain regions can be assessed using 
functional connectivity analysis, which is a method for estimating 
correlations of brain activity between regions. Resting-state 
functional connectivity analysis has been widely used to study 
the abnormalities of psychiatric diseases [18-21]. Performing 
functional connectivity analysis on resting-state data can reflect the 
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Table 1. Demographic and clinical details. 







Healthy controls 


Bipolar 1 patients 


p-value 


Number of subjects 


16 


15 




Gender (male/female) 


11/5 


10/5 


0.938^ 


Age (years)^ 


43.4±11.3 


42.6 ±9.7 


0.827'' 


HAMD^ 


1.1±1.4 


7.5±7.5 


0.003'' 


YMARSt 


0.0 ±0.0 


3.9±4.5 


0.002'' 


Duration of illness (years) ^ 




17.3±11.5 




Onset age (years) ' 




25.3 ±10.8 




Past manic episodes (times) ' 




4.3±2.4 




Past depressive episodes (times) ' 




3.0±1.4 




Drug regimen before Imaging 


Mood stabilizers 




14 




Antipsychotics 




12 




Hypnotic medications 




12 





HAMD: Hamilton Depression Rating Scale, 17 items; YMARS: Young Mania Rating Scale; 

^Continuous variables are expressed as mean ± standard deviation (SD); 

^Pearson Chi-Square test; 

''Two-tailed two-sample f-test. 

doi:l 0.1 371/journal.pone.0096422.t001 



natural mental state without producing task-induced biases and 
can facilitate investigating neural plasticity after long disease 
durations [18,19]. Previous resting-state functional connectivity 
studies on bipolar disorder have focused on anterior cingulate 
cortex connectivity [21] and the DMN [19], which are involved in 
emotional and self-referential processing. Because the functional 
connectivity regarding the FST circuit in previous studies has been 
limited by the within-circuit connectivity using emotion-relevant 
stimuli [22-24], the resting-state connectivity in the FST circuit 
and between this circuit and other regions have not been 
addressed. 

In this study, we hypothesized that the resting-state connectivity 
analysis can provide a task-unbiased observation to reveal the 
abnormahty of inherendy mental state for bipolar disorder. Based 
on the substantial studies that have reported that the bipolar 
patients exhibit deficits in mood regulations and cognitive 
functions, the altered functional connectivity may be related to 
die FST regions. Considering the facts that 1) the FST circuit was 
previously defined by the activated areas when performing 
emotional or cognitive tasks [8,9] and 2) the medial frontal 
regions of the FST circuit were primarily involved in the DMN 
during resting state [12], we anticipated that the resting-state 
connectivities related to the striatal-thalamic circuit and DMN 
may be altered in bipolar patients. To determine whether the 
striatal-thalamic circuit and DMN exist in our resting-state data, 
we used the hierarchical modular analysis to cluster the 90 resting- 
state functional regions of interest (ROIs) [25] into several modules 
according to the strength of interregional connectivity [26,27]. We 
then selected the ROIs in the modules associated with the striatal- 
thalamic circuit and DMN as the seed ROIs and subsequently 
compared the functional connectivity related to these seed ROIs 
between bipolar patients and healthy controls. 

Materials and Methods 

2.1. Participants 

The Institutional Review Board of Taipei Veterans General 
Hospital approved this study. AH participants provided written 



informed consent before participating in this study. We recruited 
15 patients with bipolar I disorder (mean age: 42.6 years; 10 men), 
characterized by at least one manic episode, from Taipei Veterans 
General Hospital. Clinical psychiatrists confirmed the diagnoses of 
all patients, using the fourth edition of the Diagnostic and Statistical 
Manual of Mental Disorders (DSM-IV). All of the bipolar patients 
exhibited no comorbidities, including schizophrenia, obsessive- 
compulsive spectrum disorders, and post-traumatic stress disorder. 
Before acquiring functional magnetic resonance images (fMRI), 
we rated the depressive and manic symptoms of bipolar patients by 
using the Hamilton Depression Rating Scale (HAM-D) [28] and 
the Young Mania Rating Scale (YMRS) [29], respectively. We 
required the bipolar patients to discontinue their medication for 1 
week before imaging. The drug regimen for the bipolar group 
before their participation in this study included mood stabilizers, 
atypical antipsychotics, and hypnotic medications. We also 
recruited 16 age- and gender-matched healthy controls (mean 
age: 43.4 years; 1 1 men) presenting no personal or family history 
of psychiatric illness. AH participants were right-handed and 
presented with no history of neurological or physical disorders, 
alcohol or drug dependence, or electroconvulsive therapy. Table 1 
lists the demographic and clinical details of all participants. 

2.2. Imaging Data Acquisition and Preprocessing 

We performed functional magnetic resonance imaging (fMRI) 
by using the Discovery MR750 3T system (GE Healthcare, USA) 
at Taipei Veterans General Hospital. We scanned each participant 
in a supine position. All of the participants were instructed to close 
their eyes, hold stiU, relax, and not fall asleep. After the resting 
scan, all of the participants reported that they could comply with 
the directions. No other functional tasks were administered to the 
participants. We acquired the resting-state fMRI of the entire 
brain in 43 axial slices with 200 volumes by using an echo planar 
imaging (EPI) sequence. The imaging parameters included a 
repetition time (TR) of 2500 ms, echo time (TE) of 30 ms, flip 
angle of 90°, field of view (FOV) of 222x222 mm, in-plane 
resolution of 3.47x3.47 mm, and slice thickness of 3.5 mm. We 
also obtained structural 3D Tl -weighted images by using a rapid 
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acquisition gradient echo for each participant (TR= 12.2 ms, 
TE = 5.2 ms, flip angle =12°, FOV = 256x256 x 168 mm^ and 
matrix size =256x256x168). 

We implemented the preprocessing of functional analysis by 
using statistical parametric mapping (SPM8; http://www.fil.ion. 
ucl.ac.uk/ spm/). We discarded the first 4 volumes of resting fMRI 
to exclude the magnetization equilibration effects and adaptation 
of the participants to the resting state. We corrected the remaining 
196 volumes for acquisition time delay among slices and realigned 
the volumes for head-motion correction (a 6-parameter alEne 
transformation). We register(-d all of tlu; corret ted images to the 
corresponding anatomical Tl -weighted images, which were 
spatially normalized to the Tl template image with a voxel size 
of 2x2x2 mm. Finally, we applied a 6-mm fuU width at half 
maximum (FWHM) Gaussian kernel to perform spatial smoothing 
on the normalized images. 

2.3. Functional Connectivity Between Pairs of 90 
Functional Regions of Interest 

We parceUated the brain into 90 functionally defined ROIs, 
based on a template proposed by Shirer et al (2012) [25]. Each 
ROI exhibited a volume of at least 30 mm'^. Using functionally 
defined ROIs is superior to using anatomically defined ROIs when 
measuring brain functional connectivity. Because the brain 
features a complex functional organization that is composed of 
small functional units, averaging time courses over functional 
parcellation is more reliable than a\eraging them over large 
anatomical parcellation, such as the Automated Anatomical 
Labeling template [25,30]. The details and abbreviations of the 
90 functional ROIs are provided in Table SI. 

We obtained the regional time series in each ROI by averaging 
the fMRI time series over all voxels in the ROI. To remove the 
bias of averaging signal to the voxels with high signal intensity, the 
fMRI time series in each voxel was normalized by its mean signal 
value. We regressed out the influence of head motion and the 
confounding effects of cerebrospinal fluid, white matter, and the 
global mean to eliminate the effects of physiological noise, such as 
respiratory-related fluctuations [31-33]. Because low-frequency 
components are the major contributors of resting-state functional 
connectivity [34], each regional time series was temporally 
bandpass filtered (0.01-0.08 Hz) to reduce the effects of low- 
frequency drift and high-frequency physiological noise. We 
calculated the functional connectivity between the filtered time 
series of each pair of ROIs by using Pearson's correlation, resulting 
in a 90 x 90 correlation matrix for each participant. We also 
applied Fisher's r-to-z transformation to ensure the normality of 
the correlation coefficients for subsequent statistical analysis [35]. 

2.4. Hierarchical Modular Analysis 

To determine whether the striatal-thalamic circuit and DMN 
exist in our resting-state data, we applied hierarchical modular 
analysis to segregate the brain regions into several modules based 
on the 90 x 90 representative correlation map. The representa- 
tive correlation map was derived from the following 3 steps: (1) the 
correlation matrices among all participants were firstly examined 
using a one-sample / test, (2) a Bonferroni-corrected significance 
level of -P<0.001 was used to threshold the resulting Pvalue matrix 
into a binary matrix for which the element was 1 if a significant 
correlation between a pair of ROIs existed and 0 otherwise, and 
(3) the values of 1 s in the binary matrix were replaced with the 
corresponding mean correlations among all correlation matrices. 

In modular analysis, modules have the characteristics that the 
connections are stronger within each module and that they are 
weaker between modules [36]. The suitability of a modular 



partition can be measured by modularity, (2,[37]: 



2m 



Im 



8{ci,cj) 



where A,j is the correlation coefficient between regions i and j, 
kj = Ay is the sum of the correlation coefficient between region 
i and its connected regions, Cj is the module to which region i is 
assigned, the 5-function d{u,v) is 1 i£ u = v and 0 otherwise (i.e., 
when regions i and j belong to the same module, the d{Ci,Cj) = 1), 

and m = ^ y . . Aij. Modularity Q, which ranges between — 1 and 

1, measures the density of links within modules compared with 
links between modules [36,37]. A larger Q, implies a superior 
partition that is more likely to constitute a modular organization 
[38].^ 

Hierarchical modular analysis is an algorithm that facilitates 
modularity optimization [26,27]. At the beginning of the 
algorithm, each region is assigned a different module. Two phases 
are then repeated. The first phase is a greedy optimization in 
which each region is designated to the same module or to an 
adjacent module. The assigned module for a particular region 
changes only if this action causes an increase in modularity. In the 
second phase, a new novel network comprising nodes that are the 
modules obtained from the first phase is constructed. A 
combination of the 2 phases is repeated until no further 
improvement of modularity is possible. When the maximal 
modularity is obtained, the optimal modular organization is 
determined. 

Based on the results of the hierarchical modular analysis, we 
selected ROIs in the modules related to the frontal, striatal, and 
thalamic regions as the seed ROIs to investigate the alterations of 
resting-state functional connectivity in bipolar patients. 

2.5. Statistical Analysis for Group Comparison 

To identify whether functional connectivity is significantly 
different from zero in each group, we applied a one-sample t test 
(P<0.001 with false discovery rate (FDR) correction for multiple 

comparisons [39]) to each functional connectivity among partic- 
ipants in the bipolar and healthy groups. For example, once 1 1 
seed ROIs related to the striatal-thalamic circuit were identified 
(see Section 3. 1 for the results of hierarchical modular analysis), 
each element of an 1 1 X 90 connectivity matrix (including the 
inter-seed-ROIs connectivit)', represented by an 1 1 x 11 matrix, 
and the connectivity between seed ROIs and other ROIs, 
represented by an 1 1 x 79 matrix) was tested by one-sample t 
test to statistically confirm the existence of functional connectivity 
for each group. After performing the one-sample t test, we 
obtained the functional connectivity exhibiting strong correlations 
(positive or negati\e; i.e., nonzero) that were consistent among 
participants in the bipolar and healthy groups. We subsequendy 
examined the significant difference of the nonzero functional 
connectivity between the bipolar and healthy groups using the 
two-sample t test as P<0.05 with an FDR correction for multiple 
comparisons [39]. 

2.6. Clinical Correlates Analysis 

We further examined the relationship between the n signifi- 
cantiy altered the functional connectivity of bipolar patients and 3 
clinical rating scales (i.e., the HAM-D, YMRS scores, and the 
duration of illness) by using a partial correlation controlling for age 
and gender effects [40]. The FDR correction for multiple 
comparisons (re X 3) was performed on the correlation between 
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Figure 1. The mean functional connectivity across all participants, (a) The nnean 90 x 90 correlation matrix across all participants. The 
anatomical locations of 90 ROIs are listed in the Table SI . The warm colors represent the positive correlations whereas the cool colors represent the 
negative correlations between ROIs. (b) The modular structure of brain functional connectivity across all participants. This modular pattern is 
obtained by reordering regions in the mean correlation matrix according to maximizing the strength of connectivity close to the main diagonal of the 
matrix. IWodule 6 consisted of 3 striatal regions and 8 thalamus subregions forming the striatal-thalamic circuit. IModule 3 consisted of 9 medial frontal 
regions, 1 temporal region, 7 parietal regions, 1 posterior cingulate cortex and 1 occipital region forming the DIVIN. The details of the modular 
structure are listed in Table S2. 
doi:l 0.1 371/journal.pone.0096422.g001 



the altered connectivity and clinical rating scales. The significant 
correlation (P<0.05; FDR correction) between each clinical rating 
scale and functional connectivity indicated that the specific 
connectivity strength might reflect the degree of depression 
(HAM-D score) or mania (YMRS score), and duration of illness 
in bipolar patients. 

Results 

3.1. Regions in the Striatal-thalamic Circuit Were 
Clustered as a Module 

Based on the results of the hierarchical modular analysis, the 90 
ROIs were segregated into 10 modules. We reordered the regions 
in the mean correlation matrix among all of the participants 
(Figure la) to illustrate the modular structure in which the strength 
of connectivity near the diagonal of the matrix was strong 
(Figure lb). Eleven ROIs, including 3 striatal regions (the ventral 
and dorsal caudate and dorsal putamen) and 8 thalamic 
subregions (the anterior, mediodorsal, ventral posterolateral, and 
pulvinar thalamus) were clustered in one module, namely, the 
striatal- thalamic circuit (Figure lb, Module 6). Nineteen ROIs, 
including 9 medial frontal regions, 1 temporal region, 7 parietal 
regions, 1 posterior cingulate cortex and 1 occipital region were 
clustered into one module forming the DMN (Figure lb, Module 
3). Because the frontal, striatal, and thalamic regions were 
separated into two modules, namely, the striatal-thalamic circuit 
and the DMN, we specified 2 sets of seed ROIs (1 1 seed ROIs for 
the striatal-thalamic circuit and 19 seed ROIs for the DMN) and 
performed one-sample t test separately. The 2 ROIs of the left 
middle temporal gyrus (Module 1) and the pons (Module 5) were 
removed in the subsequent statistical analysis because they were 
isolated from the modular structure. The details of the modular 
structure are listed in Table S2. 

3.2 Significantly Existed ROI-based Functional 

Connectivity in Bipolar Patients and Healthy Controls 

3.2.1 Striatal-thalamic-related Functional Connectivity. The 

ROI-based functional connectivities related to the striatal- 
thalamic circuit (Module 6) exhibiting significant correlations 
which were significantly diflerent from zero in each group are 
shown in Figure 2a. The first 3 seed ROIs (labeled 1, 2, and 40) 
in the y-axis of Figure 2a were the striatal regions, and the 
remaining 8 seed ROIs were the thalamic regions. Both groups 
presented significandy positive correlations within the striatal- 
thalamic circuit (Module 6). Several significandy positive 
correlations between the striatal regions and the anterior 
cingulate cortex (ROI 43 in Module 7) and between the striatal 
and thalamic regions and the posterior insula (ROI 56 in 
Module 7) were observed in both bipolar and healthy groups. 
The differences in the correlation results between the bipolar 
and healthy groups were displayed in the BD - HC difierence 
matrix in Fig. 2a. It was worthy to note that the significandy 
positive correlations in individual groups would be conservative 
and the significandy negative correlations would be artifact due 
to our preprocessing strategy using global signal regression. 



Considering the limitations of global signal regression, we only 
retained the positive correlations to construct the HCposiuve and 
BDpositive matrices. We have also displayed the differences 
between the HCpositive and BDposi,i^e niatrices (BDpojin^e " 
HCpositivc) to better illustrate the increase and decrease of 
correlations in bipolar patients without artifactual negative 
correlations (Fig. 2a). We then used the two-sample /-test as 
P<0.05 with an FDR correction for multiple comparisons to 
assess the differences of connectivities between the BDpositi,,e and 
HCposiuve matrices. 

3.2.2 DMN-related Functional Connectivity. The ROI- 
based functional connectivities related to the DMN (Module 3) 
exhibiting significant correlations in each group are shown in 
Figure 2b. The first 9 seed ROIs (labeled 6, 29, 35, 60, 21, 16, 36, 
8, and 38) in the y-axis of Figure 2b were the medial frontal 
regions, and the remaining 10 seed ROIs included the temporal, 
parietal regions, posterior cingulate cortex and occipital region. 
The significantly positive correlations within the DMN were 
observed in both groups. Both groups exhibited a similar pattern 
of connectivity between the seed ROIs in DMN and ROIs in other 
modules. For example, several significandy positive connectivities 
were observed between seed ROIs in DMN and the retrosplenial 
and posterior cingulate cortices and the precuneus (ROI 19, 15, 
and 26 in Module 4) and between seed ROIs and the middle 
frontal gyrus (ROI 84 and 80 in Module 9). The correlation 
difference matrices (BD - HC and BDposjtive - HCposMve) related to 
the DMN was also shown in Fig. 2b. 

3.3 Significant Changes in Striatal-thalamic-related 
Functional Connectivity for Bipolar Patients 

The results obtained from the two-sample / test revealed that the 
bipolar patients presented significantly decreased and increased 
connectivities (marked by white asterisks in BDp,,^;,;,,,. - HCp„s;,i„,. 
matrix in Fig. 2a) related to the striatal-thalamic circuit, compared 
with the healthy controls. The significantly decreased functional 
connectivities were between the bilateral striatum regions (ROI 1, 
2, and 40) and the middle as well as posterior cingulate cortex 
(ROI 10 and 25), and between the bilateral striatum (ROI 1, 2, 
and 40) and the mediodorsal as well as ventral posterolateral 
thalamus (ROI 57, 73, and 75). The part of anterior thalamus 
(ROI 1 and 2) also exhibited decreased connectivity with the 
middle cingulate cortex (ROI 10) as well as with the mediodorsal 
and ventral posterolateral thalamus (ROI 57, 73, and 75). 

The significandy increased functional connectivities for the 
bipolar patients were between the left thalamus subregions (ROI 
73 and 77) and the right parahippocampus (ROI 22). Table 2 
shows the results of the significandy changed connectivity in the 
bipolar patients compared with healthy controls. 

3.4. Clinical Correlates of Striatal-thalamic-related 
Functional Connectivity 

The relationships between the 9 significantly different connec- 
tivities between bipolar patients and healthy controls and 3 clinical 
measures were examined with correction for 9 x 3 comparisons. 
The results of the clinical correlates indicated that the illness 
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Figure 2. The significantly existed ROI-based functional 
connectivity in healthy subjects and bipolar patients and the 
correlation difference matrix between the bipolar and healthy 
groups. The x axis represents 88 ROls (the 2 isolated ROIs, including 
the left middle temporal gyrus and the pons were removed), and the y 
axis represent the seed ROIs. (a) Striatal-thalamic-related Functional 
Connectivity. The first 3 seed ROIs are the striatal regions and the other 
are the thalamic regions, (b) DMN-related Functional Connectivity. 
HCpositive and BDposi,i„e matrices were constructed by retaining the 
postive correlations. The significantly different connectivity in the 



HCpositive matrix was marked by a white asterisk (*). 
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duration was negatively correlated witli tlie functional connectiv- 
ities between the striatal regions and the middle cingxilate cortex, 
including tliat between the left striatal region and the middle 
cingulate cortex (R= —0.560, P= 0.047, uncorrected) and that 
between the right caudate and the middle cingulate cortex (R = — 
0.626, P= 0.022, uncorrected) (Figures 3a and 3b). Although the 
clinical correlations were insignificant after FDR corrections, the 
functional connectivities related to the striatal-thalamic circuit in 
the bipolar patients tended to decrease with the illness duration. 

3.5. Significant Changes in DiVlN-related Functional 
Connectivity for Bipolar Patients 

The results obtained from the two-sample / test revealed that no 
significant changes in frontal-related functional connectivity in the 
DMN module in the bipolar patients, compared with the healthy 
controls. Two significandy decreased functional connectivities 
(marked by white asterisks in BDp„s;ri„e - HCp„sitive matrix in 
Fig. 2b) were observed between the left middle occipital gyrus 
(ROI 18 in Module 3) and the left posterior ventral parahippo- 
campal gyurs (ROI 1 7 in Module 4) and between the left middle 
occipital gyrus and the right angular and middle occipital gyrus 
(ROI 23 in Module 3) in the bipolar patients. Based on the results 
of the one-sample / test, 201 DMN-related functional connectiv- 
ities were significantly existed in bipolar and healthy groups. The P 
value of the two-sample / test after correction for 201 comparisons 
was considerably strict, that is the significant difference was 
identified if P<0.0002 = 0.05/201. Accordingly, only 2 signifi- 
candy altered functional connectivity related to the DMN can be 
found in bipolar patients. 

Discussion 

This study investigated the resting-state functional connectivity 
changes to reveal the abnormalities of a mental state without 
producing task biases for bipolar patients. We hypothesized that 
the resting-state functional connectivity related to the striatal- 
thalamic circuit and DMN might be altered in bipolar patients. 
Based on the results of the hierarchical modular analysis that the 
striatal and thalamic regions were clustered into one module in our 
resting data, we compared the connectivity related to the striatal- 
thalamic circuit between bipolar patients and healthy controls. 
The medial frontal regions were excluded from the module of the 
striatal and thalamic regions because the medial frontal regions 
exhibited strong connectivity with several parietal, temporal 
regions and the posterior cingulate cortex to form another module 
belonging to the DMN (Module 3 in Table S2). The strong 
correlations between the medial frontal regions and the striatal 
and thalamic regions might be observed when performing 
emotional or cognitive tasks rather than during the resting state 
[8,9] . Finally, our results indicated that bipolar patients presented 
changed resting-state connectivity related to the striatal-thalamic 
circuit. The possible explanations for the altered resting-state 
connectivity of bipolar patients are discussed in the subsequent 
sections. 

4.1 Connectivity Changes in the Resting State 

Resting state data, which is measured by the absence of 
predetermined experimental task, provide several traits for 
probing psychiatric disorders. First, because bipolar disorder is 
related to internal anomalies in feeling and thinking, the resting 
state allowed bipolar patients to relax and spontaneously switch 
into self-referential representation reflecting their natural mental 
state without causing task-induced biases [18,19]. Additionally, the 
resting-state connectivity can be altered when an individual suffers 
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Table 2. Significantly changed connectivity related to the striatal-thalamic circuit in bipolar patients. 







Significantly changed functional 
connectivity 


Connectivity strength In healthy 

(mean-i-std) 


Connectivity strength In bipolar 

(mean-i-std) 


p-value (FDR corrected) 


Decrease in bipolar patients 


1 L. vCAU, dPUT, aTHAOlO 


MCC 


0.343 ±0.1 94 


0.108±0.170 


0.023 


1 L. vCAU, dPUT, aTHA O 75 


R. vpl THA 


0.341 ±0.261 


0.115±0.199 


0.049 


2 R. vCAU, dPUT, aTHA olO 


MCC 


0.306±0.172 


0.1 11 ±0.1 66 


0.031 


2 R. vCAU, dPUT aTHA o 57 


R. mdTHA 


0.926±0.304 


0.630±0.210 


0.031 


2 R. vCAU, dPUT aTHA o 73 


L. vpl THA 


0.515±0.266 


0.297±0.121 


0.038 


40R. dCAUO 10 


MCC 


0.309±0.205 


0.030±0.176 


0.013 


40 R. dCAU O 25 


MCC, PCC 


0.384±0.268 


0.068±0.229 


0.019 


Increase in bipolar patients 


73 L. vpl THA « 22 


R. pvPHIP 


0.044±0.220 


0.248±0.162 


0.042 


77 L. puTHA <:> 22 


R. pvPHIP 


0.133±0.227 


0.338±0.166 


0.038 
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from a disease with a long duration [19]. Because the bipolar 
patients in this study suffered from several episodes of abnormal 
mood disturbances for many years (average duration: 17.3 years), 
die influence of long-term mood abnormalities might alter the 
resting-state functional circuits. Therefore, we used resting-state 
data to explore the brain functional alterations in bipolar patients. 

4.2 Decreases in Striatal-related Functional Connectivity 
in Bipolar Patients 

We obsei-ved that the bipolar patients exhibited significantly 
decreased striatal-related functional connectivity, including the 
connectivity between the striatum and thalamus as well as between 
the striatum and middle and posterior cingulate cortex (Table 2). 
The abnormalities of resting-state connectivity could be associated 
with the disrupted structural fiber integrity. The striatum has 
complex fiber connections with the thalamus, which relays 
subcortical input to the cortex (i.e., the cortico-striato-paUido- 
thalamic circuit) [4 1 ,42] . Previous studies using diffusion weighted 
imaging have reported that bipolar patients exhibited a disrupted 
white matter structure (i.e., decreased fractional anisotropy to 
reflect the reduced fiber integrity) between the striatum and 
thalamus in the internal capsule [43,44]. The disruption of fiber 
connections between the striatum and thalamus in bipolar patients 
can possibly cause decreased resting-state connectivity. Consider- 
ing this basis of the structural deficits, the decreased resting-state 
functional connectivity related to the striatum and thalamus can 
be used as a probe to correlate with the task performance involving 
emotional and reward processes for bipolar patients [45]. 

Another significantly decreased functional connectivity in 
bipolar patients was observed between the striatum and the 
middle and posterior cingulate cortex (Table 2). A similar finding- 
was reported in depressive patients by using resting-state fMRI 
[46] . Another study of healthy participants determined that the 
resting-state connectivity between the striatum and the middle and 
posterior cingulate cortex was associated with dopamine concen- 
tration [47]. Dopamine dysregulatioii has been hypothesized to be 
a factor in manic and depressive symptoms in bipolar patients 
[48,49]. The dopamine dysregulation in the pathophysiology of 
bipolar disorder could be a reason for the reduced resting-state 
connectivity between the striatum and the middle and posterior 
cingulate cortex in the bipolar patients. 



4.3 Increases in Thalamus-related Functional Connectivity 
in Bipolar Patients 

Increased connectivity was observed between the pulvinar 
thalamus and the parahippocampus in the bipolar patients 
(Table 2). Both the pulvinar thalamus and parahippocampus have 
fiber connections with the amygdala [50-52]. The regions are 
involved in emotional processing, such as encoding and retrieving 
emotional memory [53-56]. Accordingly, 2 possible rationales for 
the increased connectivity between the pulvinar thalamus and the 
parahippocampus in the bipolar patients during resting state are as 
follows. First, participants were scanned at rest without external 
stimuli, and therefore, their attention was turned to the internal 
mental state involved in memory retrieval and self-perception. In 
addition, previous studies have reported that bipolar patients 
presented a tendency to easily engage with emotional memory, 
particularly negative-valence memory, rather than recalling 
memory with less emotional bias [57,58]. Second, the bipolar 
patients exliibited frequent neural overactivation of the pulvinar 
thalamus and the parahippocampus during emotional processing 
[59]. This frequent coactivations between the regions might 
increase their associations [60], suggesting an enhanced synchro- 
nization among regions in the resting state. 

4.4 Limitations and Methodological Concerns 

Several limitations of this study must be mentioned. First, 
regarding the drug regimens for the bipolar group followed before 
their participation in this study, considering ethical concerns and 
clinical status of each patient, we required bipolar patients in a 
mood stable state to discontinue their medication for only 1 week 
before imaging. A 1 -week interval from the final medication might 
not be sufficient to exclude the medication effect. However, the 
elimination half-life of the drugs for the bipolar group was below 
30 hours. The plasma concentrations of the drugs were likely 
considerably low (<2%) after 7 days, which is sufTicient to 
eliminate 97%-98% of drugs [61]. Therefore, the medication 
effects should have been minimal. Second, the correlations 
between the altered connectivity and chnical measures, including 
the illness duration and HAM-D scores, were significant only 
without correction for multiple comparisons, which might be 
attributable to the limited sample size. Further studies should 
include large sample sizes to assess the clinical correlates on 
resting-state functional connectivity in bipolar patients. Finally, 
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Figure 3. Clinical correlations between the altered functional connectivity and the clinical rating scales. 
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neuropsychological testing was not performed in this study. The 
relationships between the altered resting-state functional connec- 
tivity and cognitive and behavioral deficits of bipolar patients must 
be further explored. 

Three methodological concerns must also be considered. The 
first concern is related to the impact of head movement for each 
participant during scanning. Because a maximal translation of 
below 2 mm at each axis, and a rotation of below 2° were 
exhibited in each participant, the motion-related elfects can be 
neglected. The second concern is the elfect of global mean 
regression in the data preprocessing. Although global mean 
regression can introduce negative correlations, it can substantially 
suppress erroneously positive correlations [32]. Accordingly, we 
adopted the global mean regression in our data preprocessing and 
only retained the positive correlations to construct the connectivity 



matrices for both groups. The third concern is the a priori ROI 
selection. The conventional method of a priori ROI selection is 
subjective. To alleviate this problem, we adopted a resting-state 
fMRI template and hierarchical modular analysis so that the ROI 
selection can be more systematic. 

Conclusion 

In summary, this study determined the altered resting-state 
connectivity related to the striatal-thalamic circuit for bipolar 
patients. The significantly different functional connectivity was in 
the striatal-thalamic circuit and between the circuit and other 
brain regions consisting of the middle and posterior cingulate 
cortex as well as the parahippocampus. Changes in the resting- 
state connectivity related to the striatal-thalamic circuit might be 
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an inherent basis for the altered emotional and cognitive 
processing in bipolar patients. 

Supporting Information 

Table SI Anatomical locations, abbreviations and Brodmann 

areas of 90 functional ROIs. 

(DOC) 

Table S2 Ten modules resulted from the mean correlation 
matrix across all participants by using the hierarchical modular 
analysis. 
(DOC) 
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